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ABSTRACT 

We test estimators of the galaxy power spectrum P(k) against simulated galaxy 
catalogues constructed from N-body simulations and we derive formulae to correct for 
biases. These estimators are then applied to compute the power spectrum of galaxies 
in the Stromlo-APM redshift survey. We test whether the amplitude of P(k) depends 
on galaxy luminosity, but find no significant luminosity dependence except at abso- 
lute magnitudes brighter than Mb, = —20.3, (Hq = 100km s _1 ) where there is some 
evidence for a rise in the amplitude of P(k). By comparing the redshift space power 
spectrum of the Stromlo-APM survey with the real space power spectrum determined 
from the parent APM Galaxy Survey, we attempt to measure the distortion in the 
shape of P(k) caused by galaxy peculiar motions. We find some evidence for an effect, 
but the errors are large and do not exclude a value of (5 = Q Q - 6 /b — 1, where f2 is the 
cosmological density parameter and b is the linear biasing parameter relating galaxy 
fluctuations to those in the mass, [5p/p) gal = b(Sp/p) m . The shape of the Stromlo- 
APM power spectrum is consistent with that determined from the CfA-2 survey, but 
has a slightly higher amplitude by a factor of about 1.4 than the power spectrum of 
IRAS galaxies. 

Key words: Galaxy clustering; Cosmology; Large-scale structure 



X 



1 INTRODUCTION 

The Stromlo-APM redshift survey is a 1 in 20 sparsely 
sampled subset of 1787 galaxies selected from the APM 
Galaxy survey (Maddox et al. 1990a) to a magnitude limit 
of b] — 17.15. The survey has been used in several investiga- 
tions of the large-scale clustering of galaxies: Loveday et al. 
(1992a) have investigated galaxy counts in nearly cubical 
cells with sizes in the range 5 h~ 1 Mpc ~ 60 h~ 1 Mpc; Love- 
day et al. (1995a) have computed the two-point galaxy cor- 
relation function of the Stromlo-APM survey and Loveday 
et al. (1995b) have investigated anisotropics in the two-point 
correlation function caused by redshift-space distortions (see 
Kaiser 1987). In this paper, we present an analysis of the 
power spectrum, P(k), of galaxies in the Stromlo-APM sur- 
vey. The power spectrum of the galaxy distribution has been 
estimated from a number of other surveys (see Efstathiou 
1996 for a review). Power spectra for optically selected sam- 
ples have been computed by e.g. Baumgart and Fry (1991), 
Vogeley et al. (1992), Park et al. (1994) and Baugh and Ef- 
stathiou (1993, 1994). Power spectra for various IRAS red- 
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shift surveys have been compute by Fisher et al. (1993), 
Feldman, Kaiser and Peacock (1994, hereafter FKP) and 
Tadros and Efstathiou (1995). 

The layout of this paper is as follows. The Stromlo-APM 
redshift survey is described in Section 2. In Sections 3 and 
4, we discuss estimators of the power spectrum and their 
biases and test our techniques against simulated Stromlo- 
APM redshift surveys constructed from large N-body sim- 
ulations. In Section 5 we present estimates of the power 
spectrum for volume limited and flux limited samples se- 
lected from the Stromlo-APM survey and we investigate 
the sensitivity of the results to the volume limit and the 
weights assigned to the galaxies. A number of authors have 
claimed that luminous galaxies are more strongly clustered 
than less luminous galaxies (e.g. Hamilton 1988, Santiago 
and da Costa 1990, Iovino et al. 1993, Park et al. 1994). We 
investigate the dependence of clustering strength with lumi- 
nosity in Section 5. In Section 6 we investigate the distortion 
of the power spectrum measured in redshift space caused by 
galaxy peculiar motions and we attempt to measure the ef- 
fect by comparing the redshift-space estimates of P(k) from 
the Stromlo-APM survey with real-space estimates of P(k) 
inferred from the parent APM catalogue (see Baugh and 
Efstathiou 1993). We present our main conclusions in Sec- 
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tion 7 and compare our estimates of P(k) with those for the 
CfA-2 survey and IRAS surveys and with the predictions of 
various cold dark matter models. 



2 DATA SET 

The Stromlo-APM redshift survey is described in detail by 
Loveday et al. (1992a). The survey covers an area approx- 
imately defined by the equatorial coordinates 21 h < a < 
5 h ,-72.5° < 5 < -17.5°. In the analysis described be- 
low, we analyse flux limited and volume limited samples. 
The volume limited samples are constructed by removing all 
galaxies with redshifts z > z max and absolute magnitudes 
M > M cr u, where 



mu m — 25 — 5log (d L {z max )) — kz„ 



(1) 



mum = 17.15 is the magnitude limit of the survey 
and d,L(z max ) is the luminosity distance at redshift z max . 
Throughout this paper, we assume a spatially flat universe 
with Q. = 1, thus 



d L (z) 



2c 



(1 + *) 



(2) 



The median redshift of the Stromlo-APM survey is z = 0.05, 
hence the results presented below are insensitive to the as- 
sumed cosmological model. We adopt a k correction of 3z in 
equation (1) which is which is appropriate for the median 
morphological type in the Stromlo-APM survey in the bj 
passband (Efstathiou, Ellis and Peterson 1988). 

There are some advantages in estimating power spec- 
tra from volume limited rather than flux limited samples. 
Firstly, every galaxy in a volume limited sample carries equal 
weight, whereas weighting factors that are a function of the 
power spectrum P(k) are required to determine a minimum 
variance estimate of P(k) from a flux limited sample if the 
underlying density field is assumed to be Gaussian (see FKP 
and Section 3). Secondly, the analysis of flux limited samples 
requires a model for the mean galaxy density n(r) as a func- 
tion of radial distance. Thirdly, by estimating power spectra 
of progressively larger volume limited samples it is possible, 
given a large enough sample, to test whether the clustering 
amplitude is a function of luminosity. These points will be 
discussed in further detail in Section 4. 

The solid line in Figure ^ shows the number of galaxies 
in volume limited subsamples of the Stromlo-APM survey as 
a function of the limiting redshift z max . The total number of 
galaxies peaks at z max — 0.06, corresponding to an absolute 
magnitude limit of — 19.34^]. For most of the analysis in this 
paper we adopt a volume limit of z max — 0.06 (coordinate 
distance of x max = 172 ft~ 1 Mpc) containing 469 galaxies. 
The dashed curve in Figure |l] shows the number of galaxies 
predicted by integrating over the galaxy luminosity function: 



N (Zma X ) — <l>* V (z max ) 



(-£)<(£)« 



t Throughout this paper, we write the Hubble constant as Hq = 
lOO/i km s" 1 Mpc - 1 and quote absolute magnitudes assuming 
h = 1 unless otherwise stated. 




Figure 1. The solid line shows the number of galaxies in vol- 
ume limited subsamples of the Stromlo-APM survey as a function 
of limiting redshift z max . The dashed curve shows the expected 
distribution derived by integrating the luminosity function from 
Loveday et al. (1992a) using the Schechter function parameters 
given in the text. 



where V (z max ) is the comoving volume of the survey to 
redshift z max , and L m %n is the luminosity of a galaxy with 
absolute magnitude M CT i t (equation 1). We have used the 
luminosity function parameters from Loveday et al. (1992b) 
0* = 1.12 x 10 -2 ft 3 Mpc" 3 , a = -1.11 and M* = -19.73; 
these are the best fitting Schechter (1976) function param- 
eters for galaxies with absolute magnitudes in the range 
—22 < M < —15, uncorrected for magnitude errors. 



3 MEASUREMENT OF THE POWER 
SPECTRUM 

In this Section we establish our notation and describe our 
estimator of P(k). We imagine that the Universe is divided 
into infinitesimal cells of volume 8V containing n» galaxies 
such that the m are either or 1. The mean and variance 
of tit are therefore, 



(n») = nSV 
(nl) = (m) 



(4) 
(5) 



The observed galaxy count n (xi) at point Xj is related to 
the real galaxy count 7ij by 



n (xi) = riiW(xi) 



(6) 



where W(x.i) defines the window function of the survey. In 
general, the window function depends on the radial distribu- 
tion of galaxies in the survey, the survey boundaries on the 
celestial sphere, and the weights applied to each galaxy. In 
this section, we assume that we are analysing a volume lim- 
ited sample in which each galaxy is assigned equal weight. 
In this case, the window function W(x) is equal to unity or 
zero according to whether a particular patch of the universe 
is included or excluded from the survey. It is straightforward 
to generalize the equations in this Section to more complex 
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window functions and to include radially dependent galaxy 
weights (see FKP). 

For the Stromlo-APM survey, the window function is 
defined by an angular mask consisting of the APM Galaxy 
Survey boundaries together with rectangular holes around 
regions on the plates where there are bright stars, satellite 
trails, globular clusters, step- wedges etc (see Plate 2 of Mad- 
dox et al. 1990b). We estimate the Fourier transform of the 
window function 



WQs) = ~ I W(x)e*' x d 3 x 



(7) 



by generating a large number random points with mean den- 
sity Hr within the survey volume and computing 



W4(k) = 

n R V t-^ 



(8) 



where the sum extends over all random points located at 
positions x». The volume V is a normalizing volume which 
encloses the survey. 

We compute the Fourier transform of the observed den- 
sity field 

1 

V ^ 



w D (k) = A > nie lk ' x 



(9) 



and we define a variable with zero mean 
A(k) =ft (k)-nH r (k) ! 



(10) 



where n is the mean galaxy density. The variance of A(k) is 
related to the power spectrum P(k) of the galaxy distribu- 
tion according to 



<|A(k)| 2 ) 



V 1 

k' 



\W(k')\ 2 +y J2 I W(k-k')| 2 P(k').(ll) 



Figure |^ shows the window function of the Stromlo- 
APM survey region volume limited at z — 0.06. This figure 
shows that the window function is sharply peaked in k-space, 
falling off approximately as fc -4 . We have computed the sum 
in equation (8) by Fast Fourier transforming the random 
density field within a cubical volume of side 840 h~ 1 Mpc. We 
also plot the components of |W(k)| 2 along the principal axes 
of the cube, showing that the window function is anisotropic 
and is narrower in the a;-direction, which is aligned along the 
radial direction of the centre of the survey. 

Since the window function is sharply peaked in k-space, 
it is a good approximation to remove -P(k') from the sum- 
mation in equation (11), thus our estimate of P(k) is 



Pe(k) = 



-|A(k)| 2 -^-^ 



V 

v s 



E k ,|^e(k')l 



(12) 



where No and Nr are the total number of galaxy and ran- 
dom points within the survey volume Vs, 



Vs 



W(x) d 3 x. 



(13) 



The negative terms in equation (12) correct for Poisson shot 
noise in the galaxy and random number distributions. Notice 
that in equation (12) we have assumed that the mean galaxy 
density is equal to Nc/V s . The estimate P e will therefore be 
biased low because the estimate of the mean galaxy density 



°x n O<x 



0.01 



0.1 



k ( h Mpc 1 ) 

Figure 2. The Fourier transform of the Stromlo-APM survey 
volume limited at z = 0.06. The solid line shows the window 
function averaged in radial shells in k-space while the crosses, 
squares and triangles show \W(k x , 0, 0)| 2 , \W (0, k y , 0)| 2 , and 
|W(0, 0, k z )\ 2 components respectively. The x-axis is aligned ap- 
proximately along the radial direction pointing to the centre of 
the survey on the celestial sphere. We have normalized the curves 
so that they are equal to unity at k = 0. 



is affected by galaxy clustering. An expression for this bias 
is derived in Appendix A. 

We estimate the sampling errors on the power spectrum 
in two different ways. The first method is based on the error 
analysis of FKP who assume that the galaxy density field 
is a Gaussian point process. In the case of a volume limited 
survey, equation (2.4.6) of FKP for the error on the point 
with wave vector k is, in our notation, 



V a (l + ^+nP(k)) 



N G 



N, 



(14) 



In equation (14), k' and k" are constrained to lie in the 
same shell with modulus k and N S um is the total number of 
terms in the double summation. 

The second method is based on the variance of P(k) de- 
rived from mock Stromlo-APM surveys constructed from N- 
body simulations. The mock surveys, and a comparison with 
the FKP errors derived from equation (14) are described in 
Section 4.2. This is a useful test because equation (14) as- 
sumes Gaussian statistics, whereas the true galaxy density 
field is non-Gaussian at least on scales < 5 ft -1 Mpc where 
the rms fluctuations exceed unity. By comparing the er- 
rors from mock catalogues which are non-Gaussian on small 
scales, we can check the accuracy of equation (14) over the 
full range of wavenumbers presented in this paper. 

In the above analysis, we have assumed that we are 
analysing volume limited samples and assigning an equal 
weight to each galaxy. The analysis of a flux limited sam- 
ple is more complicated because the window function Vy(x) 
depends on radial distance, and because we must assign radi- 
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ally dependent weights to the galaxies to minimise the vari- 
ance on P(k). These points are discussed in detail by FKP 
and we follow their analysis in this paper. If the underlying 
density fluctuations are Gaussian, FKP show that the vari- 
ance on the estimated P(k) is minimised if each galaxy is 
assigned a weight 



w(r) = 



l + n{r)P(k) : 



(15) 



where n(r) is the mean galaxy density as a function of ra- 
dial distance r, which we compute from the luminosity func- 
tion of the Stromlo-APM survey with the parameters given 
in Section 2. From equation (15), we see that the mini- 
mum variance weighting for each wavenumber k, depends 
on the true value of the power spectrum at wavenumber k. 
Rather than applying different weights at each wavenumber, 
we have estimated the power spectrum from the flux-limited 
survey for four values of P(k) in equation (15), P (k) = 4000, 
8000, 16000 and 32000 ( /i _1 Mpc) 3 , which span the range of 
interest at wavenumbers < 0.3/iMpc -1 . Our analysis of flux 
limited surveys follows that of FKP except that we compute 
a (the ratio of the space densities in the real catalogue to 
that in the random catalogue, equation 2.1.3 of FKP) from 
the ratio of the sums 



E 



(1 + 4nn{n)J 3 ) 



(16) 



instead of a — Nq/Nr, where the summations are over 
all galaxies and random points. We have set 471-J3 = 
10000(ft _1 Mpc) 3 (see Tadros and Efstathiou 1995). 



4 TESTS OF THE ESTIMATORS USING 
N-BODY SIMULATIONS 

In this Section, we investigate the accuracy with which we 
can recover the power spectrum from the Stromlo-APM sur- 
vey using the methods described in the previous Section. 
There are three key aspects of the analysis that we wish to 
test: (i) the assumption that the convolution of the window 
function with the power spectrum in equation (11) can be 
replaced by a product as in equation (12); (ii) the bias in the 
estimate of P(k) caused by estimating the mean galaxy den- 
sity from the survey itself (Appendix A); (iii) the accuracy of 
the FKP error estimates (equation 14) . We test these points 
by analysing mock Stromlo-APM surveys constructed from 
N-body simulations. 



4.1 Numerical Simulations 

The numerical simulations that we use here are described by 
Croft and Efstathiou (1994). They consist of three ensem- 
bles of 10 simulations each containing 10 6 particles within a 
periodic computational box of length £ = 300 h~ 1 Mpc. The 
simulations were run with the particle-particle-particle-mesh 
(P 3 M) code described by Efstathiou et al. (1985). The simu- 
lations model gravitational clustering in a cold dark matter 
(CDM) dominated universe with scale invariant initial den- 
sity fluctuations. The three ensembles are as follows: the 
standard CDM model (Davis et al. 1985), i.e. a spatially 
flat universe with Q.q = 1 and h = 0.5 (SCDM); a spatially 



flat low density CDM universe with Qq = 0.2 and a cos- 
mological constant A = = (1 - tl ) = 0.8 (LCDM); 

a spatially flat mixed dark matter model in which CDM 
contributes Q.cdm = 0.6, baryons contribute Oj, = 0.1 and 
massive neutrinos contribute f2„ = 0.3. 

For the SCDM and LCDM simulation, the initial power 
spectrum is given by 



P(k) oc 



1+ [ak+ (6fc)5 + (ck)' 



(17) 



where v = 1.13, a = h~ 1 Mpc,b = M ft _1 Mpc and c = 
ij? h~ Mpc. Equation (17) is a good approximation to the 
linear power spectrum of scale-invariant CDM models with 
low baryon density, f2f, <C £2o (Bond and Efstathiou 1984). 
The parameter F in equation (17) is equal to tloh, thus V — 
0.5 for the SCDM ensemble and T = 0.2 for the LCDM 
ensemble. The initial conditions for our MDM simulations 
are generated from the power spectrum given by equation 1 
of Klypin et al. (1993). In the MDM models, we ignore the 
thermal motions of the neutrinos and so follow the evolution 
of a collisionless cold component with Oo = 1. 

The final output times of the models are chosen to ap- 
proximately match the microwave background anisotropies 
measured in the first year COBE maps (Smoot et al. 1992) 
ignoring any contribution from gravitational waves. Thus 
the rms mass fluctuations in spheres of radius 8 /i -1 Mpc are 
erg = 1 for the SCDM and LCDM ensembles and erg = 0.67 
for the MDM ensemble. Normalizing the anisotropies mea- 
sured from the combined first and second year COBE maps 
would increase these values of as by ~ 20-30%, but these 
small changes are unimportant for most of the discussion 
below. 



4.2 Mock Stromlo-APM Surveys 

From the simulations described above we have constructed 
mock Stromlo-APM surveys. We assume that galaxies are 
distributed like the mass points in the simulations and make 
no attempt to introduce biasing. We apply the APM angu- 
lar mask to the simulations and replicate the periodic box 
where necessary to generate distant points. We then gener- 
ate fully sampled volume limited catalogues to z max = 0.06 
that contain typically 100, 000 mass points. We also gener- 
ate volume limited samples containing on average 469 points 
by random sampling the mass points. These mock surveys 
thus have similar numbers of points as there are galaxies in 
the Stromlo-APM survey at this volume limit. Flux limited 
surveys are also generated by selecting mass points with the 
radial selection function of the Stromlo-APM survey (see 
Section 2). 

Figure |3| shows three estimates of the power spectrum in 
real space measured for each ensemble. The solid lines show 
the mean of the power spectra determined from the full cubi- 
cal N-body simulations. These are computed by Fast Fourier 
transforming the particle distribution on a 128 3 grid using 
the nearest grid point assignment scheme, as described by 
Efstathiou et al. (1985). The open squares show the power 
spectra derived from equation (12) for the fully sampled 
mock Stromlo-APM surveys volume limited at z max = 0.06. 
The symbols show the mean value of P(k) from the ten mem- 
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bers of each ensemble and the error bars show the standard 
deviation on the mean derived from the scatter in the P(k) 
estimates. The open squares are in very good agreement with 
the power spectra of the full simulations at wavenumbers 
k > 0.04/iMpc" 1 , showing that replacing the convolution 
over P(k) in equation (11) by a product is an excellent ap- 
proximation. At smaller wavenumbers, the power spectrum 
is systematically underestimated. This is caused by the bias 
associated with the mean density estimate described in Ap- 
pendix A. The dashed curves in Figure 3 show the linear 
theory power spectra for each model including an approxi- 
mate correction for the bias derived from equation (A4.2) of 
Appendix A. These curves are in excellent agreement with 
the data points, showing that the bias can be removed if the 
form of the power spectrum is known. The filled triangles 
show the power spectra for volume limited {z max — 0.06) 
realizations of the sparse sampled Stromlo-APM survey. As 
described above, these simulations have approximately the 
same number of objects as in the real catalogue and accu- 
rately model the window function of the real survey. Because 
the number of points in each realization is small (and hence 
the errors on P(k) are large), we have plotted the averages 
over 30 realizations from each ensemble. The results are very 
similar to those from the fully sampled Stromlo-APM simu- 
lations, except that the bias associated with the mean den- 
sity estimate is larger, as expected from equation (A4.2) of 
Appendix A (shown by the dashed lines passing through the 
points). 

In Figure 4 we analyse flux limited sparse-sampled 
Stromlo-APM surveys for the SCDM and LCDM ensembles 
using the techniques of FKP. Each panel shows the power 
spectra averaged over 10 realizations for four values of P(k) 
in the weighting scheme of equation (15) as indicated in 
each panel. The results are qualitatively similar to those of 
Figure 3 and show that the FKP estimator provides an un- 
biased estimate of P(k) at wavenumbers k > 0.04/iMpc _1 . 
Furthermore, Figure 4 shows that the estimates of P(k) and 
the errors on P(k) are relatively insensitive to the weighting 
scheme. This is as expected since a minimum variance esti- 
mator should not be sensitive to small departures from the 
minimum variance weighting scheme. The errors on P(k) 
at wavenumbers k > 0.1/iMpc -1 are noticeably larger if we 
adopt a constant value of P(k) > 16000(/i~ 3 Mpc) 3 in the 
weighting scheme (15), but it is clear that this is an over- 
estimate of the amplitude of the power spectrum at these 
wavenumbers. 

We have also investigated the accuracy of the error es- 
timates derived from equation (14) and its generalization 
to flux limited samples by comparing with the dispersion 
in the P(k) estimates from our mock surveys. Equation (14) 
provides an accurate estimate of the errors for wavenumbers 
k < O.l/iMpc -1 , but tends to underestimate the errors by a 
factor of ~ 2 at larger wavenumbers. 

In summary, the results of this Section show that our 
methods provide nearly unbiased estimates of the power 
spectrum at wavenumbers k > 0.04/iMpc -1 . At smaller 
wavenumbers, the estimates are biased low because we de- 
termine the mean galaxy density from the sample itself. The 
power spectra derived from flux limited samples are insensi- 
tive to the weighting scheme, provided we adopt a reasonable 
value of P(k) in equation (15). 




13.01 0.02 0.05 0.1 0.2 

k h Mpc 



Figure 3. Estimates of the power spectra in real space for the 
three ensembles of N-body simulations described in the text. The 
solid lines show the average of the power spectra derived from the 
full cubical volume of the simulations. The open squares show the 
average of P(k) derived from equation (12) applied to 10 fully 
sampled realizations of the Stromlo-APM survey, volume limited 
at z max = 0.06. The filled triangles show the average of P(k) 
over 30 realizations of the sparse-sampled Stromlo -APM survey 
volume limited at z ma x = 0.06. The error bars in each case show 
the standard deviation of the mean. The dashed lines are derived 
from equation (A4.2) of Appendix A, assuming the linear the- 
ory form of the power spectrum for each ensemble and illustrate 
how the mock surveys are affected by the bias associated with 
the estimate of the mean galaxy density. The dotted lines show 
the convolution of the linear power spectrum (equation 17) with 
the window function for the Stromlo-APM survey. The arrows 
show the smallest wavenumber represented in the simulations, 
k = 2n/ It,, where = 300 ft _1 Mpc is the computational box 
length. 



5 THE POWER SPECTRUM OF THE 
STROMLO-APM SURVEY 

5.1 Comparison of results from volume limited 
and flux limited samples 

In Figure |E| we show the power spectra for three volume 
limited subsets of the Stromlo-APM survey, together with 
la errors derived from the FKP formula (equation 14). The 
power spectra are consistent with each other, with no obvi- 
ous dependence on the volume limit. The results are quali- 
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tatively similar to those of the simulations plotted in Figure 
3. The power spectrum of the real survey has an approx- 
imately power law behaviour P(k) oc k~ 1,7 at wavenum- 
bers k > 0.05ftMpc _1 , flattens off and declines at smaller 
wavenumbers. The decline is almost certainly caused by the 
bias discussed in Appendix A and so is not a real feature of 
the galaxy distribution. The Stromlo-APM survey contains 
little information at wavenumbers k < 0.05/iMpc _1 , since 
these scales are comparable to or greater than the size of 
the survey. 

Figure 6 shows the flux limited power spectrum from 
the Stromlo-APM survey for four values of P(k) in the 
weighting function of equation (15). As in our analysis of 
the flux limited mock Stromlo-APM simulations (Figure 4) 
we see that the power spectra in Figure 6 are insensitive 
to the weighting scheme. The power spectra in Figures 5 
and 6 peak at a maximum value of ~ 30000(/i~ 3 Mpc) 3 at 
a wavenumber of ~ 0.05/iMpc _1 and decline to a value of 
~ 4000(/i" 3 Mpc) 3 at wavenumbers ~ 0.15/iMpc" 1 . Thus, 
over most of the wavenumber range plotted in Figure 6, 
a value of P(k) ~ 10 4 (ft _3 Mpc) 3 in equation (15) should 
result in an estimate of P(k) that is close to the minimum 
variance estimate. In the remainder of this paper, unless oth- 
erwise stated, we will use the flux limited power spectrum 
with P(k) set to 8000(/i _3 Mpc) 3 in equation (15). 

In Figure 7, we compare the volume limited power spec- 
trum estimate for z max — 0.06 with the flux limited estimate 
derived with P(k) = 8000(/i~ 3 Mpc) 3 in equation (15). The 
estimates are consistent with each other within the la error 
bars. This figure, and the results shown in figures 4 to 6, 
show that the power spectrum estimates from the Stromlo- 
APM survey are robust. In summary, we find no obvious 
differences between the power spectra of the volume limited 
samples plotted in figure 5, or between the power spectra 
measured from flux limited samples with different weighting 
functions. This places some constraints on variations of the 
power spectrum as a function of galaxy luminosity which we 
discuss in further detail in the next subsection. 

5.2 Variations of the clustering strength with 
galaxy luminosity 

As mentioned in the introduction, several groups have 
claimed to find various correlations between the strength 
of the clustering pattern and galaxy luminosity. In this sec- 
tion we investigate this possibility by comparing the power 
spectrum measured from volume limited subsamples of the 
Stromlo-APM data as a function of volume limit. We have 
already seen from Figure 5 that there is little evidence for 
any systematic dependence on the amplitude of P(k) with 
volume limit over the range z max — 0.05-0.07, correspond- 
ing to cuts in absolute magnitude of M crit = —18.91 at 
Zmax = 0.05 and M cr i t = —19.71 at z max = 0.07 (equation 
1). As we increase the volume limit beyond z max = 0.07, 
the number of galaxies in the Stromlo-APM survey declines 
exponentially and hence the errors on P(k) become large. 
We have therefore averaged the power spectrum estimates 
over a range of wavenumbers to test for a linear bias in P(k) 
(i.e a change in amplitude independent of wavenumber) as 
a function of volume limit. This has been done as follows: 
we have computed the power spectra of six volume limited 
samples with z max = 0.04, 0.05, 0.06, 0.07, 0.08, 0.09. These 



power spectra were averaged to produce a mean power spec- 
trum Pav(k). For each volume limit, we computed the ratio 
P(k)/P a v(k), which we averaged over the wavenumber range 
0.052 < k < 0.14/iMpc _1 to produce an estimate of the rel- 
ative bias factor b 2 (M cr u) for each volume limit. As can be 
seen from Figure 6, the error bars on P(k) depend weakly on 
wavenumber, hence it is reasonable to compute b 2 (M cri t) as 
an unweighted average over wavenumber. These bias factors 
are plotted against M cr i t in Figure 8. If there were no depen- 
dence of the clustering strength with luminosity, we would 
expect the resulting bias factors to be close to unity for all 
values of M cr u. In fact, we find that the relative bias factor is 
close to unity for most of the magnitude range, but is lower 
than unity for the faintest magnitude cut at M cr i t = —18.4 
and higher than unity for the brightest magnitude cut at 
M cr u = -20.3. 

To assess the statistical significance of these result, we 
have computed linear bias factors in exactly the same way 
for each of the 30 mock Stromlo-APM surveys constructed 
from the LCDM ensemble. The results are plotted as the 
crosses in Figure 8 together with the standard deviation for 
a single Stromlo-APM survey. The results from the mock 
surveys are close to unity over the entire magnitude range, as 
expected since there is no luminosity dependence of the clus- 
tering pattern in the simulations. These results show that 
there is no firm evidence for any luminosity dependence of 
the amplitude of the power spectrum determined from the 
Stromlo-APM survey. At the faintest magnitude cut, 7 out 
of 30 simulated catalogues gave a relative bias factor lower 
than that measured for the real data. For the brightest mag- 
nitude cut, only 2 out of 30 simulations gave a relative bias 
factor larger than that determined from the Stromlo-APM 
survey. Thus, there is tentative evidence that the amplitude 
of the power-spectrum may increase at high luminosities. 
These results seem consistent with previous work. Loveday 
el al. (1995a) found some evidence that the spatial two-point 
correlation function for galaxies in the Stromlo-APM survey 
in the magnitude range — 19 < M < — 15 has a lower ampli- 
tude than that measured for brighter galaxies. The results of 
Figure 8 suggest that this may be caused by sampling fluctu- 
ations rather than a real luminosity dependence of the clus- 
tering pattern. Park et al. (1994) have analyzed the power 
spectra of volume limited subsets of the CfA-2 redshift sur- 
vey and find evidence that the amplitude of the power spec- 
trum for galaxies brighter than M* is about 40% higher than 
the amplitude measured for fainter galaxies. This is consis- 
tent with the results plotted in Figure 8 and suggests that 
the rise in the relative bias factor at M cr u < —20.3 may be a 
real feature of the galaxy distribution. Evidently, larger red- 
shift surveys are required to firmly establish whether this 
effect is real but the results of this section show that over 
most of the luminosity range, the power spectra measured 
from the Stromlo-APM survey are consistent with the null 
hypothesis that the clustering strength is independent of lu- 
minosity. 



6 REDSHIFT-SPACE DISTORTIONS OF THE 
POWER SPECTRUM 

Galaxy peculiar velocities cause a distortion of the clustering 
pattern measured in redshift space compared to the true pat- 
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tern in real space (see e.g. Kaiser 1987). In this section, we 
analyse the distortions to the shape of the power spectrum 
measured in redshift space using the N-body simulations de- 
scribed in Section 2. This allows us to quantify the effects of 
linear and non-linear peculiar velocities and hence to estab- 
lish the range in wavenumber over which linear perturbation 
theory can be used to model the distortions. (For a similar 
analysis of redshift-space distortions in N-body simulations, 
see Gramman et al. 1993). We then compare the redshift- 
space estimates of P(k) for the Stromlo-APM survey of the 
previous section with the real-space estimates of P(k) for 
the APM Galaxy Survey derived by Baugh and Efstathiou 
(1993, hereafter BE) to quantify the effects of redshift-space 
distortion on large spatial scales. 



6.1 Redshift-space distortions of P(k) determined 
from N-body simulations 

Figure 9 shows the power-spectra measured in real-space and 
redshift space for each of the ensembles of N-body simula- 
tions described in Section 2. In each case, the amplitude of 
the power spectrum measured in redshift space is enhanced 
on large scales compared to the power spectrum measured in 
real space. In linear perturbation theory, the power spectrum 
in redshift-space, P s (k), is related to the power-spectrum in 
real-space, P r (k), according to the formula 



P s (k) = P r (k) (1 + fin 



2\2 



fi 



b ' 



(18) 



(Kaiser 1987). In this equation, p is the angle between the 
vector k and the line-of-sight, 6 is a linear bias factor relating 
fluctuations in the galaxy distribution to fluctuations in the 
mass distribution, (5p/p) ga i = b(Sp/ p) mass , and we have 
assumed that power spectrum is determined from a patch of 
the universe that subtends a small solid angle at the position 
of the observer. Averaging equation (18) over the angle p, 
we obtain 



p a (k) = p r (k) (i + |a + 



(19) 



Thus the ratio P s (k)/ P r (k) depends on the parameter fi — 
Q ' 6 /b in linear perturbation theory^]. 

On small scales (wavenumbers k > 0.1/iMpc -1 ), the 
redshift-space power spectra in Figure 9 are suppressed in 
amplitude compared to the real-space spectra. This is caused 
by the small-scale peculiar velocities in high density regions 
that produce so called 'fingers of God' in redshift space. 
The solid lines in Figure 9 show a simple model derived 
by Peacock and Dodds (1994) that incorporates both the 
linear theory distortion of equation (18) with a model for 
the distortions caused by fingers of God. The small scale 
peculiar velocities are assumed to be uncorrelated in position 
and drawn from a Gaussian distribution with 1-dimensional 
dispersion a v . The combined effects lead to the following 
formula relating the angle averaged power spectra in redshift 
space and real space: 



P B (k) = P r (k)G(fi,ka v ) 



(20) 



t This is true even for low density spatially flat models, see Pee- 
bles (1984) and Lahav et al. (1991) 



where the function G is given by 

G(fi, ka v ) = ^ [3/? 2 + 4/3 (ka v ) 2 + 4 [ko v f] 



8 {ka v f 



exp y— (fcu) 2 ) 
4 (far) 4 



2 ^3 + 2(fca) 2 ) +4/3 (kaf] 



(21) 



The solid lines in Figure 9 show equation (20) for each en- 
semble, where we have used the 1-dimensional velocity dis- 
persion determined from the simulations a v — 670.5, 344.1, 
and 510.0 km s" 1 for the SCDM, LCDM and MDM en- 
sembles respectively. These curves provide a good match 
to the results from the N-body simulations over the range 
of wavenumbers plotted in Figure 9 and are more accurate 
than the model described by Gramman et al. (1993). Nev- 
ertheless, the model of equation (20) is simplistic, as it ig- 
nores non-linear streaming of galaxy pairs (cf. eg. Nusser 
and Fisher (1995)). 

Figure 10 shows the ratio of the power spectra plotted 
in Figure 9. The solid lines show the linear theory relation 
of equation (19) and the dotted lines show equation (20). 
These results show that at wavenumbers k > 0.1/iMpc -1 , 
small-scale peculiar velocities cause a significant reduction 
in the amplitude of the power spectrum measured in red- 
shift space. One must therefore be cautious in determining 
fi from redshift distortions using the linear relations equa- 
tions (18) and (19). The range of validity of linear theory 
is strongly dependent on the amplitude of the small-scale 
peculiar velocity field, as expected from equation (20). The 
biases can be large, for example, even at k = 0.1/iMpc" 1 , 
the ratio of the power spectra in the MDM model is equal 
to 1.57 compared to the linear theory prediction of 1.87. 

The amplitude of the small-scale peculiar velocity field 
is relatively poorly constrained by observations. For exam- 
ple, Davis and Peebles (1983) found an rms relative peculiar 
velocity of V12 = 340 ± 40kms _1 between galaxy pairs sepa- 
rated by ~ 1 /i _1 Mpc from an analysis of the CfA-1 survey 
and Bean et al. 1983 found U12 = 250 ± 50kms _1 from an 
analysis of a smaller, but deeper survey. However, values as 
high as 600-800kms _1 have been derived in the literature 
(Hale-Sutton et al. 1989, Mo et al. 1993). The best con- 
straints on v\2 for optical galaxies come from an analysis of 
the CfA2 and Southern Sky Redshift Survey by Marzke et al. 
(1995) who find v i2 = 540 ± 180kms _1 . The error on V12 is 
large and the Marzke et al. results are consistent with the 
small-scale peculiar velocities (a v ~ V12/V2) for the mass 
points in the LCDM and MDM models but are lower than 
those of COBE normalized SCDM models. In principle, with 
a larger survey, we could use equations (20) and (21) (or a 
more complicated non-linear model) to simultaneously de- 
termine the parameters fi and a v (see Cole et al. 1995, for 
an application to redshift surveys of IRAS galaxies) . 

In addition to the effects of non-linear peculiar veloci- 
ties, biases in the estimates of P(k) derived from equation 
(12) will also introduce departures from the linear theory 
predictions for the ratio P 3 (k)/ P r (k) at small wavenum- 
bers. This is illustrated in Figure 11, where we have plotted 
the ratio of the redshift-space power spectra of the mock 
Stromlo-APM simulations to the real-space power spectra 
measured from the full N-body simulations. The decline in 
the ratio at k < 0.05/iMpc -1 is caused by the biases in the 
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power spectrum estimates discussed in Appendix A. Even 
at wavenumber k ~ 0.05/iMpc -1 , the ratios for the LCDM 
and MDM ensembles fail to reach the linear predictions pri- 
marily because the convolution of the redshift space power 
spectrum with the window function causes a slight depres- 
sion in its amplitude (c/ the dotted lines in Figure 3). 

6.2 Estimates of /3 from redshift space distortion 

In this Section, we estimate j3 by comparing the Stromlo- 
APM power spectrum with real space estimates of P(k) 
for the APM survey. Baugh and Efstathiou (1993) have 
described a method for recovering the three-dimensional 
power spectrum from the angular two-point correlation func- 
tion w(0). The inversion requires a knowledge of the red- 
shift distribution of the galaxies used to estimate w(9), of 
the geometry of the Universe and of the evolution of the 
power spectrum. The redshift distribution of APM galax- 
ies is well constrained by observations (see Maddox et al. 
1995). However, BE show that uncertainties in the evolution 
of P(k) and the cosmological model introduce uncertainties 
of ~ 20% in the amplitude of the real space power spec- 
trum. As an illustration of these uncertainties, we will use 
two estimates of the real space power spectrum derived from 
APM galaxies in the magnitude range 17 < bj < 20 (plotted 
in Figure 12). In these inversions we assume an Einstein- 
de Sitter universe and a power spectrum that evolves as 
P(k,z) = P(k)/(1 + z) a with a = (filled circles) and 
a — 1.3 (open circles) corresponding to a clustering pattern 
that is stable in comoving and physical coordinates respec- 
tively (see BE for further details). Uncertainties in the red- 
shift distribution and cosmological model introduce similar 
systematic differences in the real space power spectra, but 
since these errors are considerably smaller than the random 
errors in the redshift space estimates of P(k) , we will not 
discuss them in detail here. 

The open squares in Figure [l^ show the z ma x = 0.06 
redshift-space estimate of P(k) for the Stromlo-APM survey 
as calculated in Section 2.2. Qualitatively, the behaviour 
is similar to that seen in the simulations; at wavenumber 
k < 0.03 /iMpc" 1 the redshift-space power spectrum de- 
clines and falls below the real space estimates. Over the 
wavenumber range 0-03 < k < 0.07 /iMpc" 1 the redshift- 
space amplitude is enhanced suggesting a significant distor- 
tion of the clustering pattern in redshift space. At larger 
wavenumbers, the errors in the Stromlo-APM power spec- 
trum become large because it is a sparse sampled redshift 
survey and so contains little information on the small-scale 
clustering of galaxies. 

We estimate the amplitude of the redshift space distor- 
tions in the Stromlo-APM survey as follows. We compute 
the ratio P s (k) / P r (k) by linear interpolating the real-space 
power spectra plotted in Figure 12 and we average this ra- 
tio over the wavenumber range 0.05 < k < 0.1 /iMpc . 
This leads to an estimate of f3 from equation 19. To esti- 
mate the error in (3, we compute the scatter in the ratio of 
the redshift-space to real-space power-spectra from the mock 
Stromlo-APM surveys plotted in Figure 11 averaged over the 
same wavenumber range as the data. This provides an esti- 
mate of errors associated with sampling fluctuations in the 
redshift survey. We estimate the error associated with uncer- 
tainties in the real space estimate of P(k) from the scatter 





a=0.0 


a- 


=1.3 









P 


Flux Weighted P 8 (k) 


0.74 ± 0.48 


0.20 


± 0.44 


Volume Limited P s (k) 


0.38± 0.67 


-0.13 


± 0.64 



Table 1. Estimates of (3 from the Stromlo-APM survey. 

in the ratio Ps(k)/P r (k) using P r (k) determined from four 
roughly equal area zones of the APM survey (see BE). Our 
final error estimate on (3 is derived by adding these two error 
terms in quadrature on the assumption that they are uncor- 
rected. Table 1 shows our estimates of j3 derived using the 
two real-space estimates of P(k) plotted in Figure 12, and 
two estimates of the redshift-space power spectrum for the 
Stromlo-APM survey; the z max = 0.06 volume limited re- 
sults plotted in Figure 12 and the flux-limited results from 
Figure 6, with a weight function (equation 15) in which we 
have set P(k) = 8000(/i _1 Mpc) 3 . 

These results suggest a positive value of (3, though the 
errors are large. The differences in the flux-weighted and 
volume limited estimates of P(k) over the wavenumber range 
0.05 < k < 0.1 MVlpc^ 1 lead to an uncertainty of nearly a 
factor of two in the derived value of (3. A similar uncertainty 
in introduced by the difference in the real-space estimates 
of P(k) derived from the APM survey for the two adopted 
values of the evolution parameter a. 

The inferred errors on (3 are larger than those quoted 
in other analyses of redshift space distortions. For example, 
Hamilton (1993) finds (3 = 0.661^ 

from an analysis of an- 
gular moments of the two-point correlation function of the 
1.9Jy IRAS redshift survey of Strauss et al. (1992); Cole 
et al. (1995) analyse angular harmonics of the power spec- 
trum and find /3 = 0.52 ± 0.13 and (3 = 0.54 ± 0.3 for the 
1.2Jy and QDOT IRAS surveys respectively. Loveday et al. 
(1995b) find f3 = 0.48±0.12 from an analysis of anistropies in 
the redshift-space correlation function of the Stromlo-APM 
survey. At face value our estimates of (3 appear to disagree 
with those of Loveday et al. , who analyse the same survey. 
However, it seems likely that Loveday et al. have underes- 
timated the error on f3. As a further check of whether our 
errors are realistic, Table 2 list the errors and biases in (3 
derived by applying the above analysis to the mock Stromlo- 
APM surveys constructed from the simulations. Here we 
have computed (3 from the ratio of the redshift-space power 
spectra of the mock surveys to the mean real-space power 
spectrum for each ensemble, and we have neglected errors 
in the real-space estimate P(k) since these are negligible. In 
each case, the errors in f3 are close to 8/3 ~ 0.5 in good agree- 
ment with our error estimates for the real data. In addition, 
Table 2 shows that (3 is underestimated in each case be- 
cause of the nonlinear corrections described in the previous 
subsection and because of the biases in the redshift-space 
estimates of P(k) described in Appendix A. 

The results given in Table 1 are in good agreement with 
those determined by Baugh (1995) from a comparison of the 
redshift-space two-point correlation function of the Stromlo- 
APM survey and the real-space correlation function inferred 
by inverting the angular correlation function of the APM 
survey. 

Our analysis suggests that the results of Table 1 should 
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Model 





5(3 


Linear (3 


SCDM 


0.64 


0.54 


1.00 


LCDM 


0.13 


0.55 


0.38 


MDM 


0.62 


0.52 


1.00 



Table 2. Estimates of (3 derived from the ratio P s /P r over the 
wavenumber range 0.05 < k < 0.1 /iMpc -1 from mock Stromlo- 
APM surveys. 

be considered as lower bounds on the true value of (3 and 
hence that (3 = 1 is compatible with the observations. With 
larger redshift surveys it should be possible to fit a more 
complicated model to the data so extending the useable 
range of wavenumbers and reducing the biases in (3, e.g. 
equation 21 with (3 and a v as free parameters, or the non- 
linear model of Nusser and Fisher (1995) based on the Zel- 
dovich approximation. 



7 CONCLUSIONS 

The conclusions of this paper are as follows: 
[1] We have tested estimators of P{k) using simulated red- 
shift surveys constructed from N-body simulations. These 
tests show that estimates of P(k) are biased at low 
wavenumbers if the surveys are used to estimate the mean 
galaxy density. Formulae for these biases are given in Ap- 
pendix A. 

[2] We have estimated P(k) for the Stromlo-APM redshift 
survey using volume limited and flux limited samples. The 
power spectra are insensitive to the volume limit and to 
the galaxy weights applied in the analysis of flux limited 
samples. 

[3] We have investigated whether the amplitude of P(k) de- 
pends on galaxy luminosity and used N-body simulations to 
assess the statistical significance of our analysis. We find no 
evidence for any significant luminosity dependence except 
possibly at absolute magnitudes brighter than M = —20.3, 
where we find some evidence for a higher amplitude. This 
is broadly consistent with the analysis of the CfA-2 redshift 
survey by Park et al. (1994) who find that the amplitude 
of P(k) for galaxies brighter than M* is about 40% higher 
than the amplitude measured from fainter galaxies. 
[4] We have analysed numerical simulations to determine the 
effects of redshift-space distortions on the shape of P{k). 
The distortions can be well approximated by the formula of 
Peacock and Dodds (1994), equation (20), which depends on 
two parameters, (3 = O 6 /& and a measure of the small-scale 
rms peculiar velocity, a v . 

[5] We estimate the redshift-space distortions in the Stromlo- 
APM survey by comparing redshift-space power spectra 
with estimates of the real-space power spectrum determined 
by inverting the angular correlation function measured for 
the parent APM Galaxy Survey. The results indicate a pos- 
itive value of f3, consistent with other work, but the uncer- 
tainties are large and do not exclude (3 = 1. 

In Figure 13, we compare the Stromlo-APM power spec- 
trum with power spectra estimated from other surveys. The 
filled circles show the flux limited estimates of P(k) for 
the Stromlo-APM survey using P(k) = 8000 ( fr^Mpc) 3 



in the weighting scheme of equation (15). The open symbols 
show P{k) for the combined 1.2Jy and QDOT IRAS red- 
shift surveys, as analysed by Tadros and Efstathiou (1995); 
these estimates are for a flux limited sample with P(k) = 
8000 ( /i _1 Mpc)' in the weighting scheme of equation (15). 
The filled stars show P(k) from Park et al. (1994) for a vol- 
ume limited subset of the CfA-2 survey consisting of 1509 
galaxies to a coordinate distance of 101ft -1 Mpc. The power 
spectra for the two optically selected catalogues are consis- 
tent with each other within the errors, but the amplitude of 
the IRAS power spectrum is slightly lower indicating that 
IRAS galaxies are less strongly clustered than optically se- 
lected galaxies. The results shown in Figure 13 are consistent 
with a linear relative bias between optical and IRAS galax- 
ies of bopt/biRAS ~ 1-2. Recently, Landy et al. (1996) have 
computed a 2 dimensional power spectrum for the Las Cam- 
panas redshift survey and find evidence for a strong peak 
in the power spectrum at wavenumbers k ~ 0.06/iMpc _1 . 
Their results are not easily comparable to ours and it is un- 
clear whether there is any inconsistency with the three di- 
mensional power spectrum inferred for the APM Galaxy sur- 
vey (Baugh and Efstathiou 1993) or with the redshift space 
power spectrum of the Stromlo-APM survey presented here. 
A detailed comparison is deferred to a subsequent paper. 

In Figure 13b, we show the empirical estimates of P(k) 
together with redshift-space estimates of P(k) determined 
from N-body simulations of the LCDM, MDM, and SCDM 
models. The shape of the SCM curve is in good agreement 
with the power spectra measured for the IRAS and CfA-2 
surveys, though it lies low compared to the Stromlo-APM 
data points at wavenumber k ~ 0.05/iMpc _1 . The high am- 
plitude for the SCDM model implied by the COBE tem- 
perature anisotropies results in a high amplitude for the 
peculiar velocity field. This in turn leads to a significant 
difference in the shape of the redshift-space power spectrum 
compared to the real-space power spectrum (see also Bah- 
call et al. 1993) and so the model can provide an acceptable 
match to the redshift-space P{k) although it fails to match 
the shape of the real-space power spectrum (see e.g. BE, 
Maddox et al. 1995). As mentioned in Section 6.1, the small 
scale (~ l/i _1 Mpc) relative peculiar velocities predicted by 
the COBE normalised SCDM model are considerably higher 
than the relative motions between galaxy pairs measured by 
Marzke et al. (1995). Both the LCDM and MDM models 
provide acceptable fits to the power spectra measured from 
the redshift surveys and distinguishing between these mod- 
els would require accurate estimates of the effects of red- 
shift space distortions. As we have seen from the analysis 
of the Stromlo-APM survey, larger redshift surveys, such as 
the AAT 2-degree field redshift survey (see Efstathiou 1996) 
and the Sloan Digital Sky Survey (Gunn & Weinberg 1995) 
are required to determine (3 accurately and so distinguish 
between these models. 
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Appendix A: Biases in estimating P(k) 

The quantity A(fe) defined in equation (10) requires an 
estimate of the mean galaxy density n. Usually, this estimate 
will be derived from the sample itself, e.g. for a volume lim- 
ited sample we can estimate n from 

1 \ - 



Vs ^ 



niW(x.i 



(Al) 



where V s is the volume of the survey and the sum extends 
over all infinitesimal cells of galaxy count ni. Inserting (Al) 
into equation (10), we find that the expectation value of 
A(fc)| 2 is given by, 

-_> 



<|A(k)H 



i^(k')i 2 + J 



l^(k')i 2 + ^]Ti#(k-k')| 2 p(k') 



£(zi2)W(xi)W(x2)|W(k')r dVi dV 2 



■yy a I I *(zia)W r (x 1 )W(x 3 )W(k')e 



dVi dV 2 



-yy J J C(^i2)VK(x 1 )H/(x 2 )H/(k')*e- lkX2 dVi dV 2 (A2) 

The last four terms in equation (A2) account for the biases 
in estimating P(k) when the mean galaxy density is deter- 
mined from the sample itself. The first double integral in 
equation (A2) measures the excess variance in the galaxy 
fluctuations above Poisson noise averaged over the survey 
volume, 

°? = ^sj J iOnia)W(xi)W{xa) dVi dV 2 . (A3) 

If this variance is dominated by fluctuations on scales 
smaller than the scale of the survey, it is a good approxima- 
tion to set the other two integrals in (A2) equal to equation 
(A3). An approximate expression for the bias in |A(fc)| 2 is 
therefore, 



-^|W(k')| 2 -nV 2 |M/(k')| 2 

V s 



-^|H/(k')j 2 -n 2 



(A4.1) 

^|W(k')| 2 P(k') ] |W^(k')| 2 (A4.2) 
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where in (A4.2) we have written a\ in terms of the power 
spectrum P(k). We have used these equations to compute 
the dashed lines shown in Figure 3 in Section 4.2. The first 
term in each of equations (A4.1) and (A4.2) introduces a 
bias in the estimate of P(k) (equation 12) that depends on 
the number of galaxies in the survey, but is independent of 
the shape of the power spectrum; the second term introduces 
a bias that is independent of the number of galaxies in the 
sample, but depends on the shape of the power spectrum. 

It is straightforward to generalise this analysis to the 
FKP estimator of P(k) for flux limited samples. As discussed 
in Section 3.1, we estimate the mean galaxy density by com- 
puting the sum 



^— y^niw'ix^W^i) (45.1) 



Veff 

where 

Veff 



= J w'{xi)W{xi 



) dV (A5.2) 



and the weight function w'(x) is equal to 
1 



(l + 47rn(a;)J3) 



(45.3) 



Inserting this expression into equation (2.1.4) of FKP, we 
can calculated the bias in a similar way to the calculation of 
equation A2. As the final expression is lengthly, we do not 
reproduce it here. 
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Figure 4. Flux weighted power spectra determined from mock 
sparsely sampled Stromlo-APM surveys. The upper figure shows 
results for the SCDM ensemble and the lower figure shows re- 
sults for the LCDM ensemble. The solid line in each panel shows 
the mean power spectrum determined from the full N-body sim- 
ulations, as plotted in Figure 3. The points show the average 
over ten mock Stromlo-APM surveys for four values of P(k) in 
the weighting scheme of equation (15), P(k) = 4000, 8000, 16000, 
32000(/i -3 Mpc) 3 , as indicated in each panel. The error bars show 
one standard deviation on the mean. 
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Figure 5. Volume limited power spectra of the Stromlo-APM 
survey, volume limited at z = 0.05,0.06,0.07 (150 h^Mpc, 
180 /i- 1 Mpc, and 210 h^Mpc). 
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Figure 7. Comparison of the volume and flux limited power spec- 
tra estimated from the Stromlo-APM survey. The circles show the 
volume limited power spectrum for the z max = 0.06 sample as 
plotted in Figure 5 and the crosses show the flux limited power 
spectrum with P(k) = 8000(/i~ 3 Mpc) 3 in equation (15) as plot- 
ted in Figure 6. We plot la error bars, but we have suppressed 
the errors on the flux limited estimates at high wavenumbers for 
clarity. 
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Figure 6. Flux limited power spectra of the Stromlo-APM sur- 
vey for four values of P (k) in the weighting scheme of equation 
15, P(k) = 4000, 8000, 16000, and 32000(fe" 3 Mpc) 3 as indicated 
in each panel. 
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Figure 8. The relative bias factor b defined in the text plot- 
ted against limiting absolute magnitude M cr n. Circles show the 
result for the Stromlo-APM data, crosses show the results for sim- 
ulated Stromlo-APM surveys drawn from the LCDM ensemble. 
The error bars on these points show the standard deviation of one 
simulation. 
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Figure 9. Power spectra evaluated in real and redshift space 
for each ensemble. Redshift space power spectra are shown by 
the circles and real space spectra are shown by the crosses. The 
solid lines show the predicted redshift space power spectrum using 
equation 20 and the measured real space power spectrum. Values 
for cr v arc given in the text. 



Figure 10. The ratio of redshift-space to real-space power for 
the full N-body simulations (box size 300 /i^Mpc). The solid 
lines shows the linear theory predictions and the dashed lines 
show equation (20) using the same values of a v as in Figure 9. 
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Figure 12. The open squares show the redshift-space power spec- 
trum for a volume limited sample from the Stromlo-APM sur- 
vey limited at z max = 0.06. The circles show real space power 
spectra determined using the methods of BE to invert the angu- 
lar correlation function of APM galaxies in the magnitude slice 
17 < 6j < 20. The filled and open circles show the inversion 
assuming a = and a = 1.3 respectively, where the power spec- 
trum is assumed to evolve according to P(k, z) = P(k)/(1 + z) a . 
We have plotted lcr errors on the Stromlo-APM estimates, as in 
Figure 7. The error bars on the real space power spectra show the 
standard deviation determined from the scatter in the estimates 
derived from 4 roughly equal area zones of the APM survey as 
discussed in BE. 



Figure 11. The ratios of redshift-space power spectra measured 
for sparsely sampled mock Stromlo-APM surveys to the real- 
space power spectra for the full N-body simulations. The solid 
and dashed lines are as in Figure 10. 
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Figure 13. Estimates of the power spectra for various rcdshift 
surveys. The filled circles show our flux-limited estimates of P(k) 
for the Stromlo-APM survey, as plotted in Figure 7. The open 
circles show a flux limited estimate of P(k) for IRAS galaxies from 
Tadros and Efstathiou (1995). The filled stars show estimates 
of P(k) from Park et al. (1994) for a volume limited subset of 
the CfA-2 survey. The three curves in Figure (13b) show redshift 
space power spectra determined from the mass distributions in 
N-body simulations of the three CDM-like models described in 
the text; erg gives the rms amplitude of the mass fluctuations in 
spheres of radius 8/i -1 Mpc and we have multiplied the power 
spectrum of the MDM model by b 2 = (1.2) 2 to approximately 
match the power-spectra of the Stromlo-APM and CfA-2 surveys. 
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